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Abstract. It is well known that, when analyzed at the light of current synthesis model predictions, variations in the physical 
properties of single stellar populations (e.g. age, metallicity, initial mass function, element abundance ratios) may have a similar 
effect in their integrated spectral energy distributions. The confusion is even worsened when more realistic scenarios, i.e. 
composite star formation histories, are considered. This is, in fact, one of the major problems when facing the study of stellar 
populations in star clusters and galaxies. Typically, the observational efforts have been aimed to find the most appropriate 
spectroscopic indicators in order to avoid, as far as possible, degeneracies in the parameter space. However, from a practical 
point of view, the most suited observables are not, necessarily, those that provide more orthogonality in that parameter space, 
but those that give the best balance between parameter degeneracy and sensitivity to signal-to-noise ratio per A, SN(A). In 
order to achieve the minimum combined total error in the derived physical parameters, this work discusses how the functional 
dependence of typical line-strength indices and colors on SN(A) allows to define a suitability parameter which helps to obtain 
better realistic combinations of spectroscopic data. As an example, we discuss in more detail the problem of breaking the 
well known age-metallicity degeneracy in relatively old stellar populations, comparing the suitability of different spectroscopic 
diagrams for a simple stellar population of solar metallicity and 12 Gyr old. 
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1. Introduction 

1.1. Evolutionary synthesis modeling 

In order to predict the expected spectral energy distribution 
(SED) of simple stellar populations (chemically homogeneous 
and coeval stellar systems), it is possible to use first prin- 
ciples (e.g. initial mass function, star formation rate, stellar 
isochrones, element abundance ratios) to generate synthetic 
star systems. This technique, known as evolutionary synthe- 
sis modeling, has been widely employed to understand the 
origin and evolution of star clusters and galaxies (Crampin 
& Ho yle [19611 Tinsley fl972l fl978l fl980l Tin sley & Gunn 
fl976l Gunn, Stryker & Tinsley [T98T1 B ruzual 119551 120021 
Aragon-Salamanca, Gorgas & Rego 1987 Chariot & Bruzual 
[T99~n Bruzual & Charlot fT993l Arimoto & Yoshii [T986lfT987l 
Guiderdoni & Rocca-Volmerange \TW7l Buzzoni [T9^9l[T995l 
Mas-Hesse & Kunth [T99T1 Fritze-v. Alvensleben & Gerhard 
[T994I Cervi fio & Mas-Hesse fl994l Worthey fl994l W orthey 
& Ottaviani 119971 Bressan, Chiosi & Fagotto 119941 Chiosi, 
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Bressan & Fagotto 1996 Tantalo et al. 1998 , Milone, Barbuy & 
Bica [19951 Leifherer & Heckman lT995l Leitherer et al. fT^96l 
[T9^9l Fioc & Rocca-Volmerange ITW71 Vazdekis et al. [T9^6l 
uW71 I2UU51 Vazdekis [T9391 Mayya fT995l [1W71 Garcfa- 
Vargas, Molla & Bressan Molla & Garcia- Vargas [2000 
Maraston 119981 Schiavon, Barbuy & Bruzual EEM Origlia 
& Oliva 120001 Zackrisson et al. 2001] Thomas, Maraston & 
Bender 2003 1. 

The reliability of model predictions has greatly increased 
as their developers include more realistic physical ingredients. 
However, as discussed by Chariot, Worthey & Bressan 1 1996 1, 
there are still problems due to uncertainties in the theory of 
stellar evolution (e.g. post-main-sequence stages), the physics 
of stellar interiors (e.g. atomic diffusion, helium content, the 
temperature of the red giant branch), and the lack of complete 
stellar spectra libraries. It is important to note that although ini- 
tially it is straightforward to predict spectroscopic indices from 
this type of models, the inherent problems associated to the 
SED libraries, either empirical or theoretical, have a non neg- 
ligible influence in the line-strength predictions. For instance, 
empirical SED libraries constitute a coarse grained, and usu- 
ally incomplete (specially for nonsolar metallicities and non- 
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solar abundance ratios) sampling of the atmospheric stellar pa- 
rameter space, whereas theoretical libraries usually exhibit sys- 
tematic discrepancies among themselves and when compared 
with observational data (e.g. Lejeune, Cuisinier & Buser 1997, 
1998). 

The use of empirical fitting functions (e.g. Gorgas et 
al. IT9931 IT9991 Worthey et al. IT9941 Cenarro et al. 120021 
can help to reduce substantially these effects (Worthev ll994l 
Vazdekis et al. I2003i . They do not only allow the computa- 
tion of line-strength indices for any given combination of in- 
put parameters, but the error in their predictions can be mini- 
mized with the use of a large set of stars. However, and since 
the empirical fitting functions only predict the value of a given 
line-strength feature for a fixed set of stellar atmospheric pa- 
rameters, the shape of the spectrum that leads to such value 
is therefore unknown. To insert the fitting function predictions 
into the evolutionary synthesis models it is necessary to use 
the local continuum of each single star in the SED library as a 
reference continuum level. In this way it is possible to weight 
the luminosity contribution of each type of star, in the neigh- 
borhood wavelength region of each index, to obtain the final 
line-strength prediction. 

In addition, there are also additional sources of biases in 
model predictions. Cervino et al. ( I2UUU1 |2"0"0"T1 I2UU21 . and 
Cervino & Valls-Gabaud ( 2003 1 have thoroughly analyzed the 
impact of the actual discreteness of real stellar populations (see 
also Bruzual l2001l and references therein), the Poissonian dis- 
persion due to finite populations in non-time-integrated observ- 
ables, and the influence of the interpolations in time-integrated 
quantities, among others. 

But far from being a discouraging situation, the recogni- 
tion of all these problems is providing a solid understanding 
of the challenging task of modeling stellar populations. In this 
sense, the collective effort of many modelers (e.g. Leitherer et 
al. 1996 1 is given strength to the idea that reliable and unbiased 
model predictions are starting to emerge. 

1.2. Physical parameter degeneracy 

Although spectroscopic data provide a direct way to analyze 
the integrated light of composite stellar systems, the predic- 
tions from simple stellar population synthesis models reveal 
that variations in the relevant physical properties of such sys- 
tems may produce quite similar spectral energy distributions 
(SEDs). This conspiracy leads to undesirable degeneracies 
when passing from the observable space (e.g. that defined by 
line-strength indices and colors), to the parameter space (age, 
metallicity, initial mass function, etc.). 

Among the best known examples of degeneracy we must 
highlight the one exhibited by age and metallicity in the study 
of relatively old stellar populations (O'Connell I1976I I1980I 
113331 Aaronson et al. 119751 WortheyQ994 Faber et al. [T9^4l . 
This outstanding problem drove many authors to seek for 
spectral line-strength indices which were more sensitive to 
age than to metallicity and vice versa (e.g. Rose 1985 1994 
Worthey 19941- in this sense, Worthey (1994) introduced an 
interesting quantitative measure of the metal sensitivity of each 



index, computed as the partial derivatives c/log(age)/c/log(Z) 
around his model predictions for a 12 Gyr old stellar popula- 
tion with solar metallicity. Since then, large efforts have been 
focused toward the search of spectral features with very high 
(e.g. Fe4668) and very low (e.g. Hfi, Hy) metal sensitivities. 
However, this work has led to the use of individual and narrow 
absorption features (e.g. Jones & Worthey 1 19951 Worthey & 
Ottaviani 1997; Vazdekis & Arimoto 1999 1 for which accurate 
measurements demand high signal-to-noise ratios. In addition, 
these spectral signatures are usually very sensitive to spectral 
resolution and, therefore, velocity dispersion. 

1.3. Compromise between orthogonality and errors 

It is important to note that since the problem is to break a de- 
generacy, in practice the real concern is how uncertain the re- 
quested physical parameters are when derived from a particu- 
lar observable space. In this sense, two circumstances have to 
be carefully handled. The first is the orthogonality of the iso- 
parameter lines in the observable space. As we have just men- 
tioned, this is precisely the major concern of previous works. 
The second condition to be aware of is the propagation of the 
errors in the spectroscopic indices into the corresponding er- 
rors in the parameters. However, and as it is expected, narrow 
indices (better suited to provide more orthogonality) exhibit 
larger errors than broad spectral features, for a given signal- 
to-noise ratio. Summarizing, orthogonality and small errors 
are magnitudes that can not be, a priori, simultaneously maxi- 
mized. As a result, it seems clear that the most suited observ- 
able space will be that in which the two mentioned require- 
ments are best balanced. 

The relevance of finding this equilibrium can hardly be 
overemphasized, specially when one considers the important 
observational effort that is being (or is going to be) spent in 
ambitious spectroscopic surveys, like e.g. DEEP (Mould fl993l 
Ko o 119951 . EFAR (Wegneret al. Hf96t . CFRS (L illy et 
al. [19951 Ham mer et al. \TWi\ . Sloan (York et al. 120001 
Kauffmann et al. l2"U0Hl, VL T-VIRMOS (Le Fevre et al. l2"UU0l 
SDSS (Bernardi et al. 12003b . In all these type of surveys, a large 
amount of spectroscopic data is collected, although signal- 
to-noise ratios and spectral resolution are typically moderate. 
These factors strengthen the need of a quantitative estimation 
of the reliability of the physical parameters derived from such 
spectroscopic studies. 

With a clear practical sense, in this paper we explore the 
way to determine those combinations of spectroscopic observ- 
ables that provide robust tools to constrain physical properties 
in stellar populations. For this purpose, we are going to as- 
sume that evolutionary synthesis model predictions are error 
free. Although, as we have discussed in Sect. 11.11 this is not 
at present the case, we want to concentrate in the problem of 
balancing errors and degeneracy. For this reason, model uncer- 
tainties and biases are out of the scope of this paper. In sec- 
tion |2] we review and enlarge our previous work concerning 
random errors in line-strength indices, showing that a common 
functional dependence of final index errors on signal-to-noise 
ratio can be found for different index and color definitions. In 
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section |3 we obtain simple formulae to quantify total errors 
in the derived physical parameters when derived from spec- 
troscopic measurements. As an illustrative example, we exam- 
ine in more detail the age-metallicity degeneracy in section 0] 
through the comparison of the suitability of different spectro- 
scopic diagrams for a 12 Gyr old simple stellar population with 
solar metallicity. We summarize the conclusions of this work 
in section □ Finally, some more technical aspects have been 
deferred to Appendices lAHCl 



2. Errors and spectroscopic measurements 

2.1. Error handling 

Random uncertainties and biases are inherently associated 
to the physical process of data acquisition. Random errors 
can be easily derived with the help of statistical methods. 
Unfortunately, the situation is not so simple when handling sys- 
tematic effects, where a case by case solution must be sought. 
In practice, the aim is to obtain reliable quantitative constraints 
of the total random errors present in the data while having 
uncorrected systematic effects (if any) well within the range 
spanned by the former. For this to be the case, possible sources 
of systematic effects should be identified and alleviated during 
the measure process. In this paper we are assuming that this is 
actually the case, and for that reason we are exclusively focus- 
ing on the impact of random errors. 

Although, as we have just mentioned, appropriate observa- 
tional strategies can greatly help in reducing the sources of data 
biases, the unavoidable limited exposure time that can be spent 
in each target determines the maximum signal-to-noise ratio 
in practice achievable. The data reduction process, aimed to 
minimize the impact of data acquisition imperfections on the 
measurement of data properties with a scientific meaning for 
the astronomer, is typically performed by means of arithmeti- 
cal manipulations of data and calibration frames. As a result, 
the initial random errors present in the raw scientific and cali- 
bration data are combined (and thus enlarged) and propagated 
throughout the reduction procedure. 



In a recent paper, Cardiel et al. (2002 1 have discussed the 
benefits and drawbacks of different methods to quantify ran- 
dom errors in the context of data reduction pipelines. One of the 
conclusions of this work is that a parallel reduction of data and 
error frames is likely the most elegant and general approach, 
and, in some circumstances, the only way to proceed when ob- 
serving or computing time demands are specially restrictive. 
It must be noted, however, that in order to apply this method 
to compute final errors, it is essential to prevent the introduc- 
tion of correlation between neighboring pixels, which danger- 
ously leads to underestimated errors. This problem arises when 
one performs image manipulations involving rebinning or non- 
integer pixel shifts of data, which is the case of those data re- 
duction steps devoted to correct for geometric distortions, to 
produce a wavelength calibration into a linear scale, or to cor- 
rect for differential refraction variations with wavelength, to 
mention a few. Fortunately, a modification in the way typical 
reduction processes operates can help to solve this problem. 
Although we are not going to enter into details (we refer the 



interested reader to that paper), the key point is to transfer the 
responsibility of the most complex reduction steps to the analy- 
sis tools, which must manipulate data and error frames using a 
distorted system of coordinates, overriding the orthogonal co- 
ordinate system defined by the physical pixels in a detector. 

2.2. Random errors in spectroscopic measurements 

Once it can be assumed that reliable final random error esti- 
mates are available, and that, in comparison, systematic biases 
are not relevant, it is straightforward to obtain a quantitative es- 
timate of the error in a given spectroscopic measurement. Since 
the information collected by detectors is physically sampled in 
pixels, the starting point in the analysis of a single spectrum 
will be the spectrum itself SiAf) (with i = 1, -/V p i xe i s ) and its 
associated random error spectrum cr s (/I,). In the following dis- 
cussion, we are assuming 

cov[S(Ad,S(Aj)] = (S(A i )S(A j ))~(S(A i ))(S(Aj)) = Q, 

Vije [l.flpiaaJ, (1) 

i.e. errors are not correlated. From here, two different ap- 
proaches can be followed (Cardiel et al. 1998] Cardiel 1999 1. 

One possibility is to estimate numerically the effect of er- 
rors via Monte Carlo simulations. In practice, new instances 
of the spectrum, § (At), can be generated introducing Gaussian 
noise in each pixel using, for example, 



S(Ad = S(At) + V2 o- s (Ai) V-In(l -n) cos(27rr 2 ), 



(2) 



where r\ and r 2 are two random numbers in the range r\ , r 2 e 
[0, 1). The unbiased standard deviation of all the measure- 
ments, performed over a sufficiently large number of simulated 
spectra, provides the final error. 

Another method to estimate errors consists in the use of 



analytical formulae. In fact, Cardiel et al. (1998 1 and Cenarro 
et al. (1200 It have already presented analytical expressions 
to compute errors in the 4000 A break (D4000; defined by 



Bruzual 1983), and in atomic, molecular and generic indices. 
Interestingly, in the case of molecular indices, and when atomic 
and generic indices are measured in magnitudes, index errors 
can be derived by (see e.g. Appendix A in Cenarro et al. 2001 1 



cr(M) 



c(M) 
SN(A)' 



(3) 



where cr(M) is the random error in the index M, SN(A) is the 
averaged signal-to-noise ratio per A measured in the pixels in- 
cluded in all the bandpasses which define the considered index, 
and c(M) is an index dependent constant. A similar expression 
can be deduced for the D4000 (Cardiel 1999 1, namely 



c(D 40 oo) - 



D 



4000 



1 



10 SN(A) 



(4) 



However, in order to take advantage of the fact that the relative 
error in D4000 is exclusively a function of SN(A), it is more 
useful to measure this kind of spectral feature in logarithmic 
units, i.e. 



D4000 = 2.5 log 10 (D 4 ooo), 



(5) 
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since, with this definition, 



o"(D 4 ooo) 



2.5 log 10 e 1 



10 



c(D 4 ooo) 



SN(A) SN(A) 



(6) 



which has the same functional form that Eq. Q- However, 
since the definition of this index covers a relatively large wave- 
length range, the computation of the averaged signal-to-noise 
ratio in the considered bandpasses may introduce a systematic 
bias in the estimation of the random error. This point is treated 
in more detail in Appendix lAl In addition, in Appendix 151 we 
show that Eq. also holds when using colors. 

In summary, errors in typical spectroscopic measurements 
(line-strength indices and colors) can be accurately estimated 
as a constant divided by an appropriate average of the signal-to- 
noise ratio per A. TableQJlists those constants for common line- 
strength features and colors. Note that for colors we have em- 
ployed Eq. jB.16t . For classical molecular indices (and atomic 
indices measured in magnitudes), which are defined with the 
help of three bandpasses, the error coefficients are computed as 
(see Eqs. (44) and (45) in Cardiel et al. 1998 1 



c(M) = 2.5 log 10 (e)x 




(7) 



where Ah, A c and A r are the central wavelengths of the blue, 
central and red bandpasses, respectively, and AAh, AA C and AA r 
are the bandpass widths 1 . 

It is important to highlight that the coefficients c(M) ob- 
tained with the previous expression are valid as long as the 
three bandpasses do not overlap. If this is not the case, the 
coefficients can be computed numerically through numerical 
simulations. This is in fact the situation for the three narrow 
HyvA.cr indices introduced by Vazdekis & Arimoto ( 1999}, in 
which part of the central bandpass containing the spectral fea- 
ture overlaps with the blue bandpass of the pseudo-continuum. 
In Fig.H we show the result of simulating 4000 spectra using 
Eq. (E), from which we have derived the c(M) coefficients for 
HyvA.o- by fitting the measured random errors as a function of 
SN(k). 
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Table 1. Coefficients to estimate the expected random error in 
typical line-strength indices and colors, as a function of the 
mean signal-to-noise ratio per A, following Eq. 0. Note that 
these coefficients are valid for atomic indices when measured 
in magnitudes (i.e. as if they were molecular indices). For line- 
strength indices measured in A see coefficients in Table 1 in 
Cardiel et al. ([1998 1. 



prediction of evolutionary synthesis models, as a linear func- 
tion of n variations of physical parameters Pj around the point 
(p° v p\, ...,p .,..., p°„) of the form 



Am 



(8) 



7=1 



being (p^, p^, p°j, p^) the physical parameters associated 
to the observables (to", to°, ...,nfi, m"), and where in general 
the ctjj coefficients are functions of (m^mS, ...,m9, ...,m„). If 
the n x n matrix A of the coefficients a,j is not singular (i.e. 
det(A) + 0), we can also express locally variations in the phys- 
ical parameters as a function of variations in the spectroscopic 
measurements by 



3. Reliability of physical properties derived from 
spectroscopic data 

Given an n-dimensional observational space, built with the 
help of n spectroscopic measurements M,, with i - \,...,n, 
small variations in all those measurements around a given point 
{m° v m^, ...,rrfi, m°) in that space can be expressed, using the 

1 For those readers interested in computing expected errors in 
atomic indices measured in A, it is possible to use 

cr(h) - 5 — , 

SN(A) 

where c t and c 2 coefficients for typical line-strengh features are given 
in Table 1 of Cardiel et al. t'1998l. Note that in this case the error 
estimates depend on the absolute index value. 



n 

Apj a J^bji Ami. (9) 

Since each set of spectroscopic measurements will be, in 
practice, affected by random errors cr(m,) (assumed to fol- 
low a normal distribution), a random sampling of a given as- 
tronomical object provides a cloud of points of coordinates 
(m 1 ,m2, ...,m ; , m n )* (with k = 1, A^ po i nts ) scattered around 
(m?, OTj, m r m„) and following a multivariate normal dis- 
tribution. Assuming that the bandpasses defining m, and mj do 
not overlap V/, j with i + j (i.e. their random errors are uncorre- 
cted), the surfaces of constant probability density are hyperel- 
lipsoids centered at (m° v m Q 2 , rrfi, m"), with axes parallel to 
the coordinate axes. The volume Vm of the hyperellipsoids (for 
a fixed confidence level 1 - ar) in the n-dimensional space de- 
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Fig. 1. Random errors from numerical simulations in the mea- 
surement of the HyvA.o- indices defined by Vazdekis & Arimoto 
J1999i . as a function of the mean signal-to-noise ratio per A. 
The fit to the data points (full line) provides the error coeffi- 
cients c(M) defined in Eq. 0. The dashed line indicates what 
should be the location of the data if the bandpasses defining 
those indices did not overlap. 



(<*,») n 
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SN(k)i 



where 



K=|det(B)| f]e(m ; ) = 



1 



det(A)| 



i=l 



and 



(a, n) = 



2 [nxlia)] 
n T(«/2) 



n/2 



(12) 



(13) 



(14) 



Summarizing, for a fixed number of spectroscopic measure- 
ments, the <p(a, n) factor can be considered as a constant, and 
the volume interior to the surfaces of constant probability in the 
n-dimensional space constructed by the physical parameters Pj 
(j — 1, n), can be expressed as 



V P oc 



«n 



1 



(15) 



where k encapsulates both the sensitivity of the spectro- 
scopic indices to the signal-to-noise ratio — through the 
c(m,) coefficients — , and the degree of degeneracy between 
the derived physical parameters — by way of the a,- 7 or b^ 
coefficients — . For that reason k is a parameter that provides 
quantitative information concerning the suitability of a given 
set of spectroscopic indices to the study of the corresponding 
physical parameters. For our purposes, the best observational 
space will be that for which this suitability parameter attains 
the lowest value (i.e. minimum Vp volume). In this sense, an 
almost orthogonal parameter space greatly helps to reduce this 
number, although its precise value will also be strongly con- 
strained by the signal-to-noise dependence of the measured 
spectroscopic features. 



fined by the spectroscopic measurements will be proportional 
to the product of the semiaxes (Kendall ll961> 



Vm = 



2 [xxliu)]"' 2 

n F(n/2) 



\\ cr(mi), 



(10) 



where^(a) is the upper (lOOa)* percentile of a chi-square dis- 
tribution with n degrees of freedom, and F(n/2) is the complete 
gamma function. When considering the n-dimensional space 
defined by the physical parameters Pi, the above hyperellip- 
soids are transformed into new objects, which in general will 
be no longer hyperellipsoids. However, the volume contained 
within these n-dimensional objects can still be easily computed 



4. Age-Metallicity degeneracy: a case study 

As an illustrative example, we can analyze the typical prob- 
lem of breaking the age-metallicity degeneracy from two- 
dimensional diagrams built with line-strength indices and col- 
ors. 



4.1. Computing the suitability parameter 

In this case, the observational space is defined by two spec- 
troscopic measurements M\ and M%. Small variations in both 
indices around a given point (m?, m°) can be expressed as a 
linear function of age and metallicity of the form 



Vp = |det(fl)|Vtf, 

where B is the n x n matrix of the coefficients b 



(ID 



j,. In fact, re- 
placing Apj and Am, in Eq. l|9} by the errors cr(pj) and cr(m,) 
respectively, and using the result given in Eq. 0, the volume 
Vp can be obtained as 



2 [nx 2 n (a)] n/2 fi c(md 



Vp = |det(B)| 



n F(n/2) 



U SN(A)i 



AmA ^ (a u a n \( Alog(Z) 
Am 2 J U 2 i fl22/\Alog(age) 



(16) 



where, in general, the a t j coefficients are functions of and 
;«2. It is straightforward to show that 



an - S M , an, and 

021 
022 - 7 ' 



(17) 
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being S m, the metal sensitivity parameter introduced by 
Worthey (1994) for the I th line-strength index M\. It is impor- 
tant to note that the method followed by Worthey (1994, see his 
footnote number 4) is a good approximation for the derivatives 
flTog(age)/aTog(Z) as long as the approximation that the ay- 
coefficients are constant in the region where age and metallic - 
ity are being measured holds. However, since this is not always 
the case, in Appendixlclwe describe the use of bivariate poly- 
nomials to express locally the a,j coefficients as a function of 
age and metallicity, which in turn can be employed to derive 
more accurate local metal sensitivity parameters through the 
relations given in Eq. 

If the 2x2 matrix A of coefficients a, y in Eq. d!6t is invert- 
ible, one can also express locally age and metallicity variations 
as a function of variations in the line-strength indices by 



index-index diagram logM models 



f Alog(Z) 
Ulog(age) 



t>22 



Ami 

Am 2 



(18) 



where 

b u = 
bn = ■ 
b2\ = • 

b22 = 



a 22 /det(A), 
-ai 2 /det(A), 
-a 21 /det(A), 

an/det(A). 



(19) 



Assuming that the indices M\ and M 2 are independent, re- 
placing Ami in Eq. dl 81 by the random error <x(m/), and using 
the result given in Eq. the area of the 2-dimensional error- 
ellipse in the space defined by the physical parameters log(Z) 
and log(age), can be computed as 



Vlog(Z),log(age) - 

;r^( a )| d et(B)| 

1 



c(mi) c(m 2 ) 



SN(A)! SN(A) 2 
c(m\) c(m 2 ) 



(20) 



|det(A)| SN{A) X SN(A) 2 

from which it is evident that the suitability parameter is defined 

as 

1 



k = | det(B)| c(mi) c(m 2 ) 



sensitivity to 
age-metallicity 
degeneracy 



c(mi)c(m 2 ) , (21) 

sensitivity 
to SN(A) 



where the effect of the degeneracy of the iso-metallicity and 
iso-age model predictions are encapsulated in the geometric 
transformation represented by the a, ; coefficients (or by the b,j 
coefficients of the inverse transformation), whereas the sensi- 
tivity to SN(A) is controlled by the c(mf) factors. 

Thus, it is clear that the goal of achieving minimum total 
errors in both age and metallicity involve the balancing of the 
two effects. In fact, degeneracy and sensitivity to errors typi- 
cally behave in opposite sense as a function of the wavelength 
coverage of the considered spectroscopic feature. Narrow line- 
strength indices can be found to be sensitive almost only to age 
or to metallicity, whereas broad band features are less sensitive 
to noise, as can be seen in the functional dependence of c(mi) 



Hy HR vs. Fe4668 


1.3 


W94JW95 


H(5 A vs. Fe4668 


0.7 


W94,W096 


H<5 F vs. Fe4668 


0.9 


W94,W096 


Hy A vs. Fe4668 


0.7 


W94,W096 


Hy F vs. Fe4668 


0.8 


W94,W096 


HyvA,i25 vs. Fe4668 


1.5 


VA99 


HyvA,200 vs. Fe4668 


1.7 


VA99 


HyvA,275 vs. Fe4668 


1.8 


VA99 



Table 5. Hy V A,o- correspond to the [Hy + 1 /2(Fe i + Mg i)]^ in- 
dices defined by Vazdekis & Arimoto ( 1999 1. Model references 
are W94: Worthey ( fT994l : JW95: Jones & Worthey (T995i; 



W097: Worthey & Ottaviani (1997i; and V99: Vazdekis & 
Arimoto (1999 1. 



on bandpass widths (Eq. and Eq. 
indices and colors, respectively). 



IB. 161 for line-strength 

In Tables we present suitability parameters (more pre- 
cisely log [k]) computed for different combinations of common 
line-strength indices and colors. In their computation we have 
employed the predictions of simple stellar populations from 
Bruzual & Chariot (200 1) 2 . In addition, in Table[5]we have also 
determined the suitability parameter for those combinations of 
line-strength indices including narrow features, that so far have 
been considered to provide a better discrimination between 
age and metallicity: Fe4668 as metallicity indicator, and HyHR 
(Jones & Worthey [1995. . U6 A , U6 F , Hy A and Hy F (Worthey 
& Ottaviani nW7t . and Hy V A lC r (Vazdekis & Arimoto [T9"9"9T l as 
age indicators 3 . 

4.2. Two examples 

In Fig. |3 panels (a) and (c), we show two examples of typ- 
ical index-index diagrams. In both cases we have simulated 
the effect of observing a hypothetical simple stellar population 
12 Gyr old and with [Fe/H]=0.093 (corresponding to the phys- 
ical parameters of one of the points predicted by the grid of 
models), with a signal-to-noise ratio per A of 100. The simu- 
lations are displayed as clouds of points clustered around the 
model prediction for those physical parameters. Each simu- 
lated point in the index-index diagrams has been transformed 
into age and metallicity using a N - 2 bivariate polynomial 
transformation, as explained in Appendix lcl( see Fig. lC.U . As 
a guide to the reader, we give in Table|6]the numbers involved in 
the computation of the suitability parameter in the two index- 
index diagrams shown in Fig. [2] The difference between the 
physical parameters derived from each simulated point, relative 
to the values corresponding to the hypothetical stellar popula- 
tion, are represented as errors in log(age) and log(Z) in pan- 
els (b) and (d). The ellipses show the regions of 68.26, 95.44 



2 Available at 
http://www.sdss.mpg.de/sdssMPA/Spectral_Tools/ 

3 In the computation of the numbers listed in Tables |2JE1 we have 
assumed that all the spectral data are not correlated. Since some of the 
line-strength indices and colors overlap, we are implicitly assuming 
that overlapping features are obtained in independent measurements 
(i.e. different spectra and photometric colors). 
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Table 2. Logarithm of the suitability index, log[*], — see Eq. \2\\ — for the study of the age-metallicity degeneracy, computed 
for different combinations of common line-strength indices in the optical range. For the generation of these numbers we have 
employed the predictions of Bruzual & Chariot (2001) models, for single stellar populations of 12 Gyr old and solar metallicity. 
Better diagrams are those for which k are lower. In this sense, and as a guide for the eye, we have boldfaced and underlined the 
10% of these numbers with the lowest values. 
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Table 3. Same than Table[2]but for different combinations of colors. 



and 99.73% probability. The areas covered by these ellipses 
in the error space are clearly larger in the H/J-Fe4668 diagram 
than in the D4ooo-Fe4668 diagram (note that the axis scales in 
panels (b) and (d) are different). This result agrees with the 
larger (and thus worse) value of log[/<:] for the first index-index 
diagram. Note, however, that the errors are more correlated 
in the second diagram. In fact, the standard deviations around 
each physical parameter, displayed with thick error bars, are a 
relatively fair representation of the 1 <x error ellipse in panel (b), 



but not in panel (d). For that reason, the standard deviations by 
themselves are not a good parametrization of the actual uncer- 
tainty in the derived parameters. 

In practice, when one tries to answer the question of 
whether the integrated light of two stellar populations share the 
same underlying physical parameters within the error bars, the 
question translates into whether their error ellipses overlap in 
the space defined by those physical parameters. Since the prob- 
ability of overlapping decreases as the area of the error ellipses 
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Table 4. Same than Table[5]but for different combinations of line-strength indices and colors. 



becomes smaller, the suitability parameter is a direct indication 
of such probability. The presence of correlation between the er- 
rors is not critical, as far as the fake relationship, introduced by 
the presence of error correlation, is taken into account when 
studying measurements performed in different objects. In this 
sense, the use of numerical simulations may help to analyze the 
relative contribution of such error correlations to the intrinsic 
relationships between the physical parameters (e.g. Kuntschner 
et al. l200TV 

4.3. Selecting the most suitable diagram 

It is clear from the previous discussion that the best M1-M2 
diagram to disentangle physical parameters (and, in particular, 
age and metallicity) will be that for which the factor 



<Ami,M 2 - K 



1 



1 



(22) 



SN(A) l SN(A) 2 

is minimum. Although initially this can be achieved by select- 
ing the combination of spectral measurements (indices and col- 
ors) for which log[/c] is lowest, in practice one should consider 
realistic signal-to-noise ratios. For that reason, and although at 
the light of the results displayed in Tables GH5] it seems that 
color-color diagrams are the best option, in practice this is not 
necessarily the case. For instance, typical random errors in the 
measurement of colors are of the order of 0.01 mag (e.g. Babu 
& Feigelson 1996 1, which implies signal-to-noise ratios per A 
of the order of a few (see Fig. IB. 11 1. However, spectroscopic 
line-strength indices are commonly obtained with SN(A) ^10 
times larger. 

In order to compare the suitability of diagrams built with 
two line-strength indices, two colors, or one line-strength index 
and one color, let assume that 



-strength index ~ 10 SN(A) co]or . 
Under this hypothesis, Eq. ( 1221 provides 

« fin *coloij .color, 
Y color 1 .color, ~~ i Y ,indeX2 » 

^indexi ,indeX2 



(23) 



(24) 



and 

~ 10 ^'"dex.color {"><\ 
*/^index,color — l/^indexi ,index 2 • \*>3) 

^indexi ,indeX2 

Considering that the values displayed in Tables |2H5]corre- 
spond to log [k], and focusing in the best values (the ones high- 
lighted and underlined in these tables), we see that 

^color-color 



^index-index 

and 

^index-color 



a 0.01, 



— - 0.1. 



^index-index 

Finally, combining all these numbers we obtain 

■ best 



(26) 



(27) 



.best 

r color, .color- 



,best 

t indexi ,color-> 



,index9 " 

(28) 

This result, derived after assuming the crude hypothesis stated 
in Eq. J23b . indicates that the best index-index, index-color and 
color-color diagrams are, initially, almost equally suited for the 
study of the age-metallicity degeneracy (in 12 Gyr-old simple 
stellar populations with solar metallicity). In a real situation, 
the selection of line-strength features, colors, or a combination 
of both to build a good diagram will depend on the available 
signal-to-noise ratio for each spectral indicator. 

Focusing on the results displayed in Tables and leav- 
ing in a second plane the relevance of the signal-to-noise ratio 
just discussed, the D4000, Fe4668 and Mg2 features are the best 
line-strength features to be included in index-index and index- 
color diagrams, whereas for color-color diagrams the lowest 
log[/f] values are obtained for colors involving well separated 
bandpasses, like (g - K s ), and (V - K). 

Interestingly, the suitability parameter for combinations 
of narrow line-strength features with Fe4668 (displayed in 
Table |5Jl are worse than the value for the D4 oo-Fe4468 dia- 
gram, and only H6a, H5p, Hy A , and Hy F can rival with Mg2. 
If, in addition, we consider that the c(M) coefficients (Tabled) 
for the narrow indices are larger than the same coefficients for 
the broader features D4000 and Mg2, at a fixed signal-to-noise 
ratio per A the diagrams of these two latest spectral features 
with Fe4668 provide more information. 
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Fig. 2. Panel (a): Hy8-Fe4668 diagram, with Bruzual & Chariot (2001) models over-plotted (dashed lines and dotted lines cor- 
respond to the predictions for constant age and metallicity, respectively; ages are labeled in Gyr, and metallicities are given as 
[Fe/H]). Dots correspond to 2000 simulations of a 12 Gyr old object with solar metallicity, assuming SN(A) = 100. Panel (b): 
Errors in log(age) and log(Z) for the same simulated data displayed in panel (a). Ellipses indicate the regions of 68.26, 95.44 and 
99.73% probability, whereas the central error bars show the unbiased standard deviation in each axis. Panels (c) and (d): Same 
than panels (a) and (b) after introducing the D4000 index instead of H/3, and using the same SN(A) in the simulations. The D4000 
index is plotted in an inverse scale to display the model predictions with roughly the same orientation in age and metallicity than 
in panel (a). Note also that the axis scales in panels (b) and (d) are different. See discussion in Sect. 14.21 



4.4. Some words of caution 

It is very important to keep in mind that, for several reasons, 
the above results must be taken with care: 

(i) The list of line-strength indices and colors explored is, 
obviously, not complete. 



(ii) The derived numbers do not take into account the un- 
certainties in the stellar population modeling. In fact, in Tabled 
we compare the suitability parameters derived using the pre- 
dictions of three different sets of models. The differences ex- 
hibited by some combinations of line-strength indices reflect 
the existing discrepancies between different models. However, 
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0.069191 
0.060252 

0.051545 


0.060212 
0.052139 
0.044401 


Fe4668 


8 Gyr 
12 Gyr 
20 Gyr 


0.039794 
0.043234 
0.045004 


0.071434 
0.076271 

0.081293 


0.118949 
0.129336 
0.142160 



diagram Eq. 
parameter H/J-Fe4668 E>4ooo-Fe4668 reference 



Pio 


-0.430095 


-0.430095 


O 


Pu 


0.051432 


0.051432 




Pto 


0.040343 


0.040343 


<C2l 


Poi 


0.323917 


0.323917 


<C2l 


Phi 


-0.014719 


-0.014719 




qw 


-0.058764 


0.488931 


<Q2l 


?n 


0.003947 


-0.013260 


<Q2l 


qio 


0.002449 


-0.013883 


OJ 


qm 


-0.717729 


4.921280 


<C2l 


qoi 


0.033402 


-0.225252 


<C2l 


log[age] 


10.079182 


10.079182 




log[Z] 


0.093200 


0.093200 


<Q2l 


a n 


0.095814 


0.095814 


<16>. <C.2t 


"12 


0.031999 


0.031999 


fT6l.<C2l 


"21 


-0.018528 


0.352696 


fT6l.<C2l 


fl 2 2 


-0.044036 


0.379328 


fT6l.<C2l 


c(mi) 


0.2757 


0.1086 


EJ 


c(m 2 ) 


0.2235 


0.2235 


ED 


K 


16.995 


0.9685 


ED 



Table 6. Numerical values for the example of Section l4~2"l ('see 
also Fig0. The upper table lists the input line-strength in- 
dices employed in the computation of the bivariate polynomial 
transformations in Figs. |2ja) and|2c) around the model pre- 
dictions for a SSP of 12 Gyr with solar metallicity (boldfaced). 
The lower table displays the different parameters involved in 
the calculation of the previous transformations, and the factors 
leading to the suitability index in Eq. J21> . 



it is worth noting that for those combinations involving Fe4668 
the agreement is very reasonable, and less good for D4000, Mg2, 
andNa5895. 

(iii) The integrated light of the stellar systems under study 
(e.g. star clusters, galaxies) are not necessarily well described 
by simple stellar populations. New models, including more 
complex star formation histories, may provide a different suit- 
ability ranking of line-strength features and colors. Anyway, 
the procedure here described is still valid as long as the correct 
models are employed. 

(iv) When comparing several objects in a diagram, their 
physical parameters are expected to exhibit a range. Since the 
k values shown in Tables |2jj5] have been computed for a fixed 
age and metallicity, the suitability parameter will be different 
for each object. For illustration, in Fig.[3]we represent a three- 



2.0 




Fig. 3. Three-dimensional representation of the log[/c] values 
for all the intersecting points in the grid of models from Bruzual 
& Chariot (2001), for the H/3-Fe4668 diagram. 

dimensional plot with the log [a - ] values derived for each inter- 
secting point constituting the grid of models in the Hj0-Fe4668 
diagram. It is clear from this figure that the suitability parame- 
ter is a function of the location in the index-index plane. 

(v) Although systematic effects have not been considered in 
the analysis performed so far (see Section l2~TV in the real world 
they can be the ones that, if not properly constrained, may de- 
termine the actual suitability of a given set of spectroscopic 
observables. In practice the nature and relevance of systematic 
errors affecting line-strength indices (e.g. spectral resolution, 
radial velocity, flux calibration) and colors (e.g. zeropoint ca- 
libration, ^-corrections) are typically different. For that reason, 
and assuming that the spectroscopic data (either spectra or col- 
ors, but not both) are obtained under homogeneous conditions, 
it is safer to use the numbers displayed in Tables [2 and [3]in a 
differential way, i.e. as a guide to compare different diagrams 
within each of these tables, but not to compare the absolute val- 
ues of the suitability indices between the two tables. The same 
reasoning makes the numbers displayed in Table@]more uncer- 
tain that the ones shown in Tables |5]and[3] 

Summarizing, and even considering all these problems, it 
is clear that once a given set of models has been adopted, the 
factor displayed in Eq. \22\ is an excellent tool to estimate an 
optimized combination of spectroscopic measurements, in or- 
der to face the study of integrated spectra. 

5. Conclusions 

We have investigated the combined role of parameter degener- 
acy and signal-to-noise ratio in the study of the integrated spec- 
troscopic properties of astronomical objects. In particular, we 
have examined the effect of random errors at the light of stellar 
population model predictions. We have shown that the expected 
random error in the measurement of line-strength indices and 
colors can be very easily computed as a constant divided by an 
appropriate average of the signal-to-noise ratio per A, as stated 
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index 



models 



D4000 W94 
D4000 V00 
D4000 BC01 



Fe4668 W94 
Fe4668 V00 
Fe4668 BC01 



Mg, W94 
Mg : V00 
Mg 2 BC01 



Na5895 W94 
Na5895 V00 
Na5895 BC01 



2 

U 

0.2 
0.3 
0.2 



Z 

o 



o 



o 



0.2 1.7 0.5 
0.4 0.9 1.3 
0.4 1.2 1.9 



0.4 
0.4 
0.3 



0.7 
0.7 
0.6 



0.8 
1.1 
0.8 



0.0 
0.0 
0.0 



0.3 
0.4 
0.5 



0.4 
0.5 
0.4 



0.3 
0.7 
0.2 



0.2 
0.4 
0.0 



1.0 0.6 
0.7 0.8 
0.5 0.6 



0.4 
0.6 
0.4 



0.5 
0.7 
0.5 



0.5 
0.6 
0.6 



0.6 
0.8 
0.7 



0.4 1.0 0.6 
0.6 0.8 1.1 
0.3 1.4 0.9 



0.0 1.1 1.2 1.2 1.1 1.3 1.6 1.7 ... 1.2 2.2 1.0 0.8 1.2 1.7 1.9 1.9 2.4 3.0 1.3 1.6 1.4 
0.0 0.9 0.9 1.1 1.3 1.2 1.7 1.5 ... 1.2 2.1 0.8 0.7 1.3 1.5 1.7 1.7 3.1 2.2 1.1 1.6 1.3 
0.0 0.9 0.9 1.2 1.3 1.4 1.7 1.6 ... 1.2 2.4 0.9 0.8 1.3 1.6 1.9 1.8 2.3 3.1 1.2 1.5 1.2 



0.2 1.5 1.4 1.4 
0.4 1.2 1.4 1.2 
0.0 1.9 1.4 1.3 



1.0 
1.5 
1.3 



1.5 
1.3 
1.4 



1.9 
1.5 
1.7 



2.2 
2.4 
2.5 



0.8 
0.7 
0.8 



0.9 
0.9 
1.0 



1.3 
1.2 
1.3 



2.0 
1.7 
1.8 



1.6 1.6 
1.8 1.7 
2.8 1.8 



1.3 1.4 1.3 
1.3 1.5 1.3 
1.3 1.5 1.3 



1.4 
1.5 
1.5 



1.5 1.5 2.0 
1.7 1.4 1.4 

1.6 1.5 1.4 



0.4 
0.6 
0.3 



2.1 2.6 1.6 
2.0 2.6 1.5 
1.6 1.5 1.5 



1.3 
1.8 
1.6 



3.1 
1.9 
2.1 



2.9 
2.0 
2.5 



2.6 
2.8 
2.6 



1.3 
1.1 
1.2 



1.3 
1.3 
1.3 



1.8 
1.7 
1.7 



1.7 
1.6 
1.6 



1.5 
1.7 
1.6 



1.7 

3.7 
2.0 



2.4 
2.5 
2.7 



2.0 
1.8 
1.8 



2.1 
2.0 
2.0 



1.8 
1.8 
1.7 



1.9 
1.9 
1.9 



1.8 2.2 

1.7 1.7 

1.8 1.6 



Table 7. Comparison of log[/e] for the study of the age-metallicity degeneracy, computed from the prediction of three sets of 
models, W94 (Worthey 1994), V00 (Vazdekis et al. [T#6l Blakeslee et al. l200fl . and BC01 (Bruzual & Chariot 2001). In all the 
cases we have derived the suitability parameters for SSP around 12 Gyr old and solar metallicity. 



in Eq. Q. This simple expression allows to define a suitability 
parameter which combines both effects (degeneracy and sen- 
sitivity to noise), providing an immediate tool to compare the 
usefulness of different observational diagrams. The recipe to 
perform such a comparison is the following: 

(i) Chose a reliable stellar population model. 

(ii) Obtain the geometric transformation to convert 
the observational parameters (line-strength indices and col- 
ors) into physical data (age, metallicity, initial mass func- 
tion,. . . ). For this purpose, bivariate polynomial transforma- 
tions (Appendix ICt are a very convenient way to obtain the 
coefficients given in Eq. ( IC.2L or those corresponding to the 
inverse transformation, Eq. dl9> . 

(iii) Use Eq. J2 It to obtain k. The coefficients which indi- 
cate the sensitivity of the spectroscopic data to the signal-to- 
noise ratio can be extracted from Tabled 

(iv) Using reasonable estimates of the expected signal-to- 
noise ratio per A in each spectroscopic measurement, employ 
Eq. J22I to obtain if/. As discussed in Sect. 14. 31 this is the factor 
to be minimized. 

Classical atomic line-strength indices must be measured as 
molecular indices in order to apply the above procedure. The 
same also holds for generic indices (e.g. CaT, PaT, CaT*; see 
Cenarro et al Uffift . 

We have illustrated this method by studying in more detail 
the well known age-metallicity degeneracy. Using model pre- 
dictions for a 12 Gyr old simple stellar population with solar 
metallicity, we have shown that a broad spectral feature like 
the D4000 can be as well suited (or even more) than Hfi to an- 
alyze this kind of degeneracy, once the dependence on signal- 
to-noise ratio is taken into account. 

For all the reasons mentioned in Sect. 14.41 the aim of this 
paper is not to give a definite answer to the question of which is 
the best observational space to disentangle physical properties 
of stellar populations, but to provide an easy way to determine 
the relative suitability of different spectroscopic diagrams to 
obtain physical information of the astronomical objects under 
study. 
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Appendix A: Estimating random errors in the D 40 oo 

In order to explore in more detail the validity of Eq. @, we 
have compared the error estimations in the D4000 derived from 
that formulae with the results derived by using a more accurate 
method. 

In Fig. IA.1I we represent the relative errors in the D4000, 
e r [D4oooL as a function of the signal-to-noise per A, as mea- 
sured in the 713 spectra (including repeated observations) of 
the stellar library gathered by Gorgas et al. {1999 1 to derive 
the empirical calibration of this spectral feature. In panel (a) 
we present the measured relative error determined using an 
accurate method — Eqs. (38)-(40) of Cardiel et al. dTWgl — . 
for four different estimations of the signal-to-noise ratio per 
A. In particular, crosses and small dots indicate the results ob- 
tained when using the mean SN(A) in the blue and red band- 
passes, respectively; the signal-to-noise ratio for the open cir- 
cles has been computed as the arithmetic mean of the two pre- 
vious values, whereas for the filled circles we have employed 
the weighted mean (l/<SV>bi ue + l/(SN) red )/(l/(SN} 2 blue + 
1/(X/V)jr ed ). The prediction of Eq. @ is the diagonal full line, 
whereas the residuals with respect to this prediction are plotted 
in panel (b). The employed stellar sample typically contains 
spectra with poorer signal-to-noise ratio in the blue bandpass 
of the D4000 than in the red bandpass, and the simple arithmetic 
mean of the SN(A) is not a good approximation. 

It is clear from the previous figure that the weighted mean is 
the best approximation. We have also checked that a weighted 
mean of the form [2/(l/<SJV>£ lue + l/(SN}l d )] 1/2 (not shown in 
the figure) gives acceptable results. 

Appendix B: Estimating random errors in colors 

With the aim of estimating the dependence of random errors 
in colors on the signal-to-noise ratio, we are following here a 
similar procedure to that employed in Cardiel et al. ( 199]S|l to 
derive the corresponding formulae for classical indices. 

Given two filters, a spectral energy distribution f(A), and 
the SED of a reference object T{A) (e.g. the SED of orLyr 
for magnitudes measured in the Vega system, T{v) = 3.63 x 
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Fig.A.l. Relative errors in the D4000, ft[D400o]> as a function 
of the signal-to-noise per A, measured in the stellar library em- 
ployed by Gorgas et al. ( 1999) to derive the empirical calibra- 
tion of this spectral feature. See text in AppendixlAlfor details. 



10 20 erg cm 2 s 1 Hz 1 for AB magnitudes, or T{A) = 3.63 x 
10 9 erg ctrT 2 s _I A -1 for HST magnitudes), a color can be de- 



termined by (see e.g. Fukugita et al. 1 19951 
jf(A)R 1 (A) dA 



C = -2.5 logi 



+2.5 log 10 



Jf(A) R 2 (A) dA 
fr(X)Ri(X) dA 



I 



T(A) R 2 (A) dA 



(B.l) 



where Ri(A) is the response function of the / th filter. The previ- 
ous expression can be rewritten as 



C = -2.5 lot 



to 



(f(A)h 

<fWh 



+ C , 



(B.2) 



where 



(fW)t 



L f(A) Ri 



(A) dA 



(A) dA 



C = -2.5 log 10 — + Co 



k 2 



k 
and 



(A) dA 



(B.4) 



(B.5) 



Co = +2.5 log 



JrWRiW dA 



T(A) R 2 (A) dA 



(B.6) 



In practice integrals must be replaced by summations of the 
form 



f(A) RiA) dA^& V f{A Uk ) RU Uk ) 



(B.7) 



where is the linear dispersion (in A pixel -1 , assumed to be 
constant along the spectrum), Nj is the number of pixels cover- 
ing the z™ filter, Ays is the central wavelength of the k th pixel of 
the / th filter, and /(/lyO is the flux is that pixel. 

If (Tf(A) is the error spectrum associated to f{A), and if we 
assume that T{A) and R/{A) are error free, the expected error in 
the color can be expressed as 



2.5 log 10 e ({f(A)h 

C(C) = C 



if (A))! 

(f(A)h 



(f(A)h 



(B.8) 



The square of the last term in the previous expression can be 
written as 



Mf(A)h 

^ 2 «/(A)>i) 



((f(A)h) 2 <r 2 «/(A)> 2 ) 



(B.9) 



((f(A)h) 2 ((f(A)h) 2 ((f(A)h) 2 ' 
In the other hand, from Eq. (IB.2> it is immediate to show that 

((f(A)h) 2 



((f(A)h) : 



(B.10) 



The first fraction in the right hand side of Eq. iB.9\ can also be 
rewritten as 



= «/W)i) 2 a^OfWh) = 
((f(A)h) 2 ((f(A)h) 2 ((/(A))!) 2 

-o.4 (C -c„)l 2 ^ 2 «/W)i) 
J ((f(A))i) 2 ' 



(B.ll) 



= [10- 



It is not difficult to show that assuming that the error spectrum 
is roughly constant within each bandpass filter, making use of 
Eqs. JB.3I > and JB.5i . and after replacing the integrals by sum- 
mations, 



(B.3) cr 2 «/0l)>,-) 



«f(A))d 



<f(A)y kJ 



6- 



1 



ft- 



(B.12) 
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Fig. B.l. Random errors from numerical simulations in the 
measurement of three colors in the 131 stellar spectra from the 
library of Pickles (1998 1, as a function of the signal-to-noise 
ratio per A. The full lines are the predictions of Eq. ( IB. 16i . See 
text for details. 



where (f{A)) ^ and (cr. s ) are the mean flux and mean error 

per A in the / th filter, respectively, SN(A)i the mean signal-to- 
noise ratio per A in the i lh filter, and 

Ni 



k=\ 



k=\ 



(B.13) 



Introducing the result of Eqs. dB.9> -( fB"T2l into Eq. ( IB.8I . we 

obtain 



cr(C) = 2.5 log 10 e 
ft 



ft 



[SN(A)\] 2 [SN(A) 2 ]< 



1/2 



(B.14) 



If 5W(A)i * SN(A) 2 , this last expression adopts the same form 
that Eq. Q, 



cr(C)- 



SN{A) 
where 

c(C) = 2.5 log 10 e x (ft +ft)'/ 2 
- 1.086 (ft +ft) 1/2 . 



(B.15) 



(B.16) 



Numerical values of for typical photometric bands are given 
in TablelBTl 

For illustration, we compare in Fig. IB. li the predictions of 
Eq. ( IB.16I >. for three sample colors, with the results of numer- 
ical simulations. For this purpose, we have employed the 131 
stellar spectra from the library of Pickles ( 1998), which con- 
tains SEDs with an ample range of spectral types and luminos- 
ity classes. For each of these spectra, we have built a synthetic 
error spectrum by randomly choosing a given SN(A). Color er- 
rors were then measured in simulated instances of the spectra 



generated with Eq. 0. The full lines in Fig. IB. II are not fits 
to the data points, but the predictions of Eq. 



( IB. 16l > using the 
corresponding ^, parameters (extracted from Table lB.lY 

Finally, it is also possible to express the coefficients as 
a function of the filter width. In fact Eq. (IB. 1 3i is the discrete 
expression of the more general definition 



6 = 



f 

J filter. 



R(A) 2 dA 



j 

J filter. 



R(A) dA 



(BA7) 



If now we assume that the filter transmission can be approxi- 
mately described by a box function of the form 



R(A) = 



Rq (constant) for Ai < A < A 2 , 

for A < A\ or A > A 2 , 



it is immediate to show that 



f = (A 2 -Ai)- 1 =FW(A)-\ 



(B.18) 



(B.19) 



where FW(A) is the filter width in A. In practice, since filter 
response functions are not exactly box functions, filter widths 
can be approximately computed as 



FW(A) ( R(A) dA, 

J filter, 



(B.20) 



being R{A) the response function normalized to unity. To check 
the validity of these approximations, in Fig. IB.2l we graphically 
compare the measured filter widths versus the £ values derived 
using Eq. ( IB. 1 31 . for a collection of 360 filter response func- 
tions collected from the literature (see table caption). The full 
line is a least squares fit, whereas the dashed line is the pre- 
diction of Eq. ( IB. 191 . As it is apparent from the figure, there 
is an excellent correlation between the £ coefficients and the 
filter width, although the linear fit to the data points (full line) 
indicates that, on average, 



i * (0.78 ± 0.09) FW(ky 



(B.21) 



where the error in the coefficient is the residual standard devia- 
tion around the fit. 



Appendix C: Measuring age and metallicity from 
index-index diagrams 

Iso-metallicity and iso-age lines in index-index, index-color 
and color-color diagrams are usually far from displaying very 
regular grids, but they typically exhibit unevenly spaced and 
distorted patterns. For this reason, the computation of ages and 
metallicities from a given pair of spectroscopic measurements 
should be addressed through the use of local mapping func- 
tions which properly accounts for the geometric distortions. 
Obviously, the same is also true for the computation of the local 
derivatives (i.e. metal sensitivity parameters). 
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Table B.l. Numerical values of the ^ coefficients computed from Eq. (IB.13> . for a set of common photometric bands correspond- 
ing to the filters given in Table 9 of Fukugita et al. d!995> . Filter widths were determined with Eq. (IB.20> . 



An excellent approach to this problem is the use of bivari- 
ate polynomial transformations of the form (see e.g. Wolberg 
1992) 

N N-i 

m = X Z p v ^s^Y (iog[z]y, 

;=o ;=o 

(C.l) 

N N-i 
i=0 7=0 

where ,/V is the polynomial degree, and py and qij are the poly- 
nomial coefficients. The inverse transformation can be written 
by an analogous expression. 

For N - I the above equations only account for affine trans- 
formations (i.e. translation, rotation, scale and shear). This lin- 
ear approximation is valid when the 6 polynomial coefficients 
are derived from control points (those for which line-strength 
indices, ages and metallicities are given by the models) which 
are very close to the point (poo, goo)- Although most evolution- 
ary synthesis models provide these close control points when 
considering the line-strength predictions as a function of age, 



the same is not true for the indices variations as a function of 
metallicity. This problem leads to systematic uncertainties in 
the index predictions, as shown in Fig. lC.ll for the HyS-Fe4668 
diagram. 

The second-degree approximation, N - 2, improves the 
quality of the prediction allowing for a selection of more distant 
control points. In this case, 12 coefficients must be computed 
by solving two systems of 6 linear equations. It is straightfor- 
ward to show that the coefficients of the A matrix in Eq. dl6> 
can be rewritten as a function of the bivariate polynomial coef- 
ficients as 

an =Pio+Pu log[age] + 2 p 2 o\og[Z], 
an = Poi + P\\ log[Z] + 2 p 01 log[age], 
an = qw + qu log[age] + 2 q 2 o\og[Z\, 
an = <?oi + qu log[Z] + 2 q {)2 log[age]. 

In Fig. IC 1 1 we also show the mapping obtained for N - 2 
(thick solid lines), which clearly improves the result obtained 
for N - 1. Higher order polynomials (yV = 3, ...) are typically 
unnecessary, since the required number of points to perform the 
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Filter width (A) 

Fig. B.2. Comparison of the ^ values as a function of the filter 
width, computed with Eqs. JB. 1 31 and ( IB.20L respectively. The 
measurements have been performed using the response func- 
tion of a set of 360 photometric bands collected from the lit- 
erature (including those which are incorporated in the source 
code of the GISSEL models, and HST, KPNO and ESO filters, 
among others). The solid line is Eq. (IB.2H . a least square fit 
to the data, and the dashed line is the prediction of Eq. dB.19> . 
The difference between both lines is due to the fact that filter 
response curves are not exactly box functions. 

fit increases rapidly, and in these situations it is always possible 
to constraint the fit to a smaller region with N — 2. 

Obviously, the same procedure is valid when reading other 
physical parameters from other grids predicted by stellar pop- 
ulation models. 
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Fig. CI. Example of fit of local polynomials to describe the 
local behavior of line-strength indices as a function of physi- 
cal parameters. Dashed and dotted lines are the predictions of 
Bruzual & Chariot (2001) models for single stellar populations 
of fixed age and metallicity, respectively (ages are given in Gyr, 
and metallicities as [Fe/H]). The thin solid lines indicates the 
resulting fit after using Eq. ( IC. II with N - I, around the model 
predictions for a SSP of 12 Gyr and solar metallicity (the co- 
efficients were determined using a least-squares fit to 5 points: 
the grid model point chosen as the origin of the local transfor- 
mation, and the two closest points in both age and metallicity). 
It is clear that the fit can not be extrapolated very far from the 
central point. The thick solid lines show the result for N - 2 
(derived from the fit to the 9 points: the 5 points previously em- 
ployed in the fit for N — I, plus the 4 additional corners of the 
region delineated by the thick line). In this case it is clear that 
the N - 2 polynomial approximation provides a very reason- 
able representation of the geometric distortions when moving 
from the observational to the physical parameter space. 



